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Abstract . The objective of this paper Is to utilize the design methods of 
modern control theory to realize a "dual-adaptive" feedback control unit for 
a highly non-linear single spool airbreathing turbojet engine. 

Using a very detailed and accurate simulation of the non-linear engine at the 
data source, linear operating point models of unspecified dimension are Iden- 
tified. Feedback control laws are designed at each operating point for a 
prespecified set of sampling rates using sampled-data output regulator theory. 

The control system sampling rate is determined by an adaptive sampling algo- 
rithm in correspondence with turbojet engine performance. The result is a 
"dual-adaptive" control law that is functionally dependent upon the sampling 
rate selected and environmental operating conditions. Simulation transients 
demonstrate the utility of the dual-adaptive design to improve on-board com- 
puter utilization while maintaining acceptable levels of engine performance. 


I- INTRODUCTION 

To generate required propulsive thrust levels 
and maintain structural and thermodynamic 
engine integrity, control system designs for 
gas turbine engines of the past several decades 
only required the proper mar,.' pul at ion Of engine 
fuel flow (and possibly far. nozzle area). 
Control objectives emphasized component reli- 
ability and design simplicity to maintain 
operating performancB specifications. However, 
recent aircraft operational requirements such 
es those imposed for supersonic flight, VTOL 
capability and improved fuel conservation dic- 
tate the development of gas turbine engines 
with increased performance capabilities over 
an expanded operating environment. Conse- 
quently, turbofan engine control designs to 
meet these performance demand schedules will 
require significant increases in control flex- 
ibility. This increase in flexibility may 
assume the form of extensive variable geome- 
tries. variable pitch fan blades and/or the 
utilization of digital congiuters for engine 
performance evaluation and control. 

In this paper, we assxsne the pretence of an 
on-board digital computer to perform system 
management tasks which include the monitoring 
and contrcl of turbofan engine performance. 

In this regard, the computer may not be totally 
dedicated to engine control except perhaps 
during critical operating periods. Thus '^he 
data monitoring and control aampling rates may 
varv in correspondence with the operating 


conditions end performance requirements of the 
overall aircraft system. 

From an engine control viewpoint, the sampling 
rate and feedback gains are integrally related 
if operating point engine performar.ee is to be 
maintained at acceptable levels. To allow for 
this degree of flexibility in engine control, 
linear models (at each operating point) for a 
preselected eet of sampling rates need to be 
identified. Once the models are known, feed- 
back control laws can be obtained using sampled- 
data output regulator theory. The engine con- 
trol system sampling rate at a particular 
operating point is then determined by an adap- 
tive sampling algorithm in correspondence with 
turbojet engine performance- Repeating this 
procedure for several operating conditions 
results in a "dual-adaptive* control law that 
is functionally dependent upon the sampling 
rate selected and the environmental operating 
conditions. 

In tha next section, the Tse and Weinart [1] 
method is used to identify sodel order and 
Structure for eoveral operating point models 
of a turbojet engine over a finite aet of 
system sampling rates. The models are used in 
Section III to develop a dual adaptive control 
configuration- In Section IV the results of a 
dual adaptive control design for a turbojet 
engine are discussed. 


II. ENGINE MODEL IDENTIFICATION 






Pr-vious work in identifying gas turbofan 

er.gir.e has beer, dcr.e by Otte and Taylor 

[2], Crooks and Willshire [3], Ahibeck [4), 
Cher. [51 , Michael and Farrar [6 ] and Mueller 
[7]. Tr.ese techniques assure a model order 
to obtain linear models fror transfer function 
data [4,i,-; or by least square curve fitting 
procedures It]. In [6] the least squares 
ret'r.t d ir. (€1 was combined with a dynamic non- 
linear filter tc identify gas turbine dynamics 
frtr stochastic input output data. Using a 
perturbational rethod with an output error 
evaluation procedure, DeHoff ar.d Kali [9] 
octair. a*, iterative equation fror a synthe- 
sized tine derivative of the state vector to 
identify avatar order [10]. Although these 
nethoda are s.itable for hydror.echsr.ical con- 
trol designs, rest are not extendable to the 
digital control format established in the 
previous section. Fcr this purpose the Tse 
and Weinert nethod is selected. 
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the index > is then increntnted until d^CyM-O. 
Set p « j*-l which represents the smallest 
non-negative integer such that row i of the 
observability matrix ( A . C) is linearly depend- 
ent on the vectors in all preceeding rows. 

Since 

n 

l P. - n (5) 

i-1 


Using output data, the Tse anc Weir.ert method 
identifies a constant, multivariable, stochas- 
tic linear system which has unknown dimension, 
matrix parameters ar.d r.cise covariances. 

Since the output data set determine an equiva- 
lence class of systems [12, a general stochas- 
tic system model cannot be identified from 
steady state output data alone. However, each 
system model in this equivalence class has a 
steady state Kalman filter with a common 
impulse response and innovations covariance 
[11]. If the system matrices are chosen in a 
certain unique conor.ical form., consistant 
Kalmar, filter parameter estimates car. be 
obtained. Thus througn the identification of 
the Kalmar, filter model for the equivalence 
class, a state model for the unknown system 
car. be obtained. 


Using the notation in [1], the steady-state 
Kalmar, filter for the equivalence class is 
given by 


Vi ■ A *k * K “k 

y k - c *k * \ 


(1) 


where z^ is an n vector, y is an m vector, 

the Gaussian innovations process u^ has unknown 
covariance 


where m is the number of output variables, the 
unknown system order, n, can be determined. 
Following [1, 12] the elements in A and C can 
be determined from. (2), (3) and p^ in (5] [1]. 

This idantifieation method was applied to 
stochastic data obtained fror. an accurate non- 
linear digital simulation of a turbojet engine 
[13]. The simulation incorporates a real aas 
combustion model and experimentally determined 
compressor and turbine data. Dynamics ate 
represented by formal one-dimensional inviscid 
continuity, momentum and energy approximations 
to unsteady compressor, combustor and exhaust 
nozzle conditions. The simulation represents 
a single rotor turbojet with each steady-state 
operation point uniquely determined by rotor 
speed, E, flight math number, M, and flight 
altitude, n. 

Eight engine variables (m-B) were selected as 
system outputs: rotor speed (S) , compressor 

discharge temperature (T^) and pressure (P^), 

nozzle inlet temperature (T^) and pressure (P^) , 

turbine inlet temperature (T^) and pressure 
(F^l , and engine thrust IF). For this study 
the engine input is fuel flow, u?. 

The engine is modeled as 


and K is the Kalmar, gain matrix. Tha unknown 
vector 6. ■ {A,C,K,g} ir. addition to the model 
order n is identifiable If (1) A is stable, 
(2) the set (A,K,C) is a minimal realization 
and (3) the specific cononical structure fcr 
A (as given in [1]) is used. 

With V N » * s t * le observed data set, 

time series analysis can be used to approxi- 
mate the output covariance R(c) using 

i N ; e t 

R(C) - - > y(i*3> y (i) (2) 

N i-1 

Forming the identifiability matrix 


x - f (x, uj { > S) * ; (t) 
y ( t) - g(x, w f , s) ♦ n(t) 


( 6 ) 


n 8 

where x e R and y c R • The sampled-data 
linearized operating point siodel for (6) is 
given by 


6 Vi - A 6% ♦ » «*f * : k 
6y k " c 6 *k 4 D H 4 \ 


(7) 


where A, B, C, D, model order n and noise 
covariances are unknown. 


Data for the identification method was obtained 
from the non-linear simulation over a fifteen 
second interval (after steady state conditions 
were achieved) by perturbing the system via 
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with white noiie- Designating the engine 
rotor apeed of 36. 960 rpm as the 100% oper- 
ating point, four operating points correspond- 
ing to 60%. 90%, 100% and io4.5% rotor speed 
were selected. The 60% ar.d 104.5% operating 
points correspond to minimum and maximum rotor 
design speeds, respectively. Typical data 
histories for S and T^ are shown in Fig. 1. 

Following the procedure outlined above, R(a) 
was calculated at each operating point using 
the appropriate data set. The ident in- 
ability matrix was formed from Rio) with 
d i '(‘r) in (4) subsequently evaluated. Figure 
2 indicates d.J(yj as y varies when j»l 
(rctor speed) and j«6 (turbine temperature). 
Using a structure criterion c * 0.1. it 
fellows from Fig. 2 that py ■ 2. Evaluating 

(3) ar.d (4) for i-2,3 8 yields {p. } B . - 

l 1-2 

{0} and suggests that each output for i > 2 
is a linear combination of the two state 
variables associated with p.. Thus from (5), 
n - 2. i 

The system matrices A and C were subsequently 
evaluated at each operating point and appear 
in Table 1. To complete the model in (7), B 
ar.d D were identified using a gradient search 
procedure. With n ■ 2 and m « 8, the ter. 
unknown elements were obtained using a sum of 
squared errors cost function for gradient 
evaluation over a five second transient. The 
test input in each case was 
3 

5_,(k) - A T Sin (y k ♦ tf, ) (8) 

q*l n ’ 

with • 1 , V • 0 » Uj«.2, u.-2 and 

u. _**20 • The input magnitude, A, was adjusted 
to limit perturbations to the linear region 
about each operating point. Table 2 indicates 
the matrices B and D obtained for each operat- 
ing point. 


as a function of a system parameter to improve 
computer utilisation and sampling efficiency. 
Here sampling efficiency is defined as the 
ratio of system performance to the number of 
samples required to achieve that performance, 
l.e. , 


n 


A 


J A N A 


(9) 


where is the value of the cost function 
(inverse of performance) and N is the number 
of samples. As decreases for constant 

r\ increases, thereby improving computer utili- 
sation. 


Using (9) as a general guide, the adaptive 
sampling law [141 used in the turbojet appli- 
cation [12) is 


where 
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T . and T are lower and upper bounds on 
min max 

the sampling rate respectively. The variable 
e^ was selected as 




k 




( 12 ) 


where is engine core speed. Core speed was 
selected since it is readily measurable and a 
good indicator of dynamic engine performance. 
Sampling limits were established as 


The identification results were verified using 
the fuel flow test input of Fig. 3. The 
input is a combination step, ramp and parabola 
large enough to drive the engine operating 
point from 85% to 100% spied and was applied 
co the non-linear simulation and the compo- 
site model of Fig. 4. The composite model 
was generated by scheduling the identified 
A, B, C and D matrices as a function of rotor 
speed, S. Linear interpolation was used to 
obtain dynamics between Identified operating 
point models. Trajectory’ comparisons for S 
and T t are given in Figs. 5 and 6. Discrep- 
ancies in the figures are attributed to the 
high frequency content in T t and T (not 
shown) and to the limited spectrum r of the 
model (approximately 5 Hz.) 

In the next section, the identified operating 
point models are used to obtain constant out- 
put feedback gain matrices at each operating 
point. 


III. DUAL ADAPTIVE CONTROL 
Adaptive sailing varies the sampling race T 


T , ■ .001 sec. 

min 

T ■ .025 sec- 
max 

The variable a in (11) reflects the relative 
weighting between sampling cost (efficiency) 
and the squared error cost due to sampling. 

For the application in this paper a ■ 4(10 4 ). 

For a constant sampling rate the corresponding 
sampled-data system and steady-state feedback 
gain matrix obtained from discrete output 
feedback regulator theory [12] is time- invari- 
ant. With adaptive sampling, the corresponding 
sampled-data system is functionally dependent 
on T^ and therefore will be time-varying even 
thouah the original continuous process is time- 
invariant. If T^ were known for all k over 
the interval of control, the time evolution 
of the discrete system would be known and the 
feedback gain matrix could be obtained from 
the time-varying regulator equations . How- 
ever, {T^} is not known a priori but is a 
function of a system parameter and will vary 
with engine performance and environment- Thus 
the feedback gain matrix cannot be determined 


ir. the usual wav since at future sampling 
instants is not available. To overcome tnis 
difficulty, it is assured that T), determined 
at instant X will be sustained for all future 
tire instants. Kith this assumption, a feed- 
back control matrix car. be computed using 
steady state discrete optimal feedback regu- 
lator theory. Nov, feedback gains can be 
computed off-line for a sufficient nx-ber of 
sampling periods to establish a functional 
relationship between system gams ar.i T . 

This relationship can easily be stored in a 
digital computer memory and used on-line to 
generate the control input to the continuous 
process. 

As the system progresses through the operating 
range, the functional gain relationship be- 
comes dependent upon the operating point 
conditions as well as the sampling interval. 
Thus the dual adaptive nature cf the control 
is created. 


XV. APPLICATION [12] 

In this section the dual adaptive control of 
Section III is applied to the operating point 
models obtained in Section II. 

Case 1 y(t) ■ S 

The control task Is to accelerate the engine 
from steady-state at 90% to steady-state at 
100% rotor speed. The sampling rate was held 
constant at T ff i n with regulator weighting 
matrices selected as Q<*1C ar.d R»l. The 
resultant cost function J is then used to 
evaluate sampling efficiency t»1.'J K t and 
provides a baseline comparison for^the re- 
maining cases (set Table 1). 

Case 2 ylt) ■ S 

Adaptive sampling was incorporated into the 
control configuration of case 1. The results 
show approximately a 10% increase in J . and a 
100% improvement in sampling afficiency with 
no visible differences in the time plots. 

Case 3 y (t) ■ S 

Tc examine the effect of a change in the 
weighting coefficients of case 2, set £*50 
and retain R«l. The results are plotted with 
those of case 2 in Pigs. 7 - 10. Case 3 
control, as expected, accelerates the engine 
faster than the case 2 control. Rotor speed 
is within 1% of steady state in 0.6 sec. and 
the thrust acceleration rise time is 0.4 sec. 
as compared to 1.1 sec. and 1.2 sec. for case 
2 respectively. The penalty for this ingjroved 
acceleration of case 3 is hiqher turbine inlet 
temperature and additional fuel flow require- 
ments. 

Case 4 


the enqlne from 90% to 104.5% design speed 
without violating the 1900* R constraint. 

This temperature constrained acceleration ic 
accomplished by initially controlling the 
engine as in case 2 with y(t! • S until the 
1913‘R limit is exceeded. Then the control 
configuration is changed to include both rotor 
speed, F, and T T feedback, i.*., y(t) ■ 

|S with 0"diag [.1, .0001] and R«l. 

Once the temperature constraint is satisfied, 
the control reverts to the original case 2 
feedback arrangement. The results are docu- 
mented in Table 3 and are plotted with those 
of case 2 in Figs. 11 - 13. 


V. CONCLUSIONS 

System model identification, discrete optimal 
output regulator feedback and adaptive sampling 
methods were combined to develop a dual adap- 
tive control technique. This technique waa 
applied to a non-linear turbojet angina to 
obtain a feedback control over multi-operating 
point conditions. The example application 
clearly demonstrates the utility of the design 
approach. 
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TA5LE 1 Identified Values of A and C vs- \ Speed-Discrete Model 


» Speed Matrix Matrix Elements 
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.0127 

.0019 

-.0336 

.00029 

-.028 

.0016 

.018 
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-.00046 

-.00009 

.00418 

-.00002 
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-.00011 

-.0012 
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-.340 
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r T 
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.0144 

.0023 

-.0264 

.00089 

-.023 

.0021 

.023 
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- . 00069 

-.00017 

.oo3eo 
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.0037 

-.00017 

-.0023 
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-.256 
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1.060 








_7 
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.0153 

.0025 

-.0200 

.00051 

-.017 

.0023 

.026 
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-.00064 

- . 0 32! 

.00366 

-.00006 

.0038 

-.00023 

-.0034 

104.5 
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0 

-.318 

1 

1.119 
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.0163 

.0027 

-.018 

.00053 

-.016 

.0025 

.028 
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-.0013 

-.00049 

.0047 

-.00008 

.0036 

-.00034 

-.0044 


TABLE 2 Identified Values of B end D vs. % Speed 
» Speed Matrix Matrix Elements 


s: 

b t 

48004.7 

27815.9 







d t 

-12227.78 

-15.25 

12.47 

3508.81 

2.598 

3991.66 

13.96 

-11703.98 

90 


43653.5 

24165.8 







d t 

-10573.77 

-33.38 

4.550 

2920.50 

2.970 

3427.02 

9.0856 

114.41 

100 

b t 

45947.6 

19977.63 








d t 

-18929.90 

-223.59 

-26.322 

2535.71 

-.9140 

2870.89 

-23. 63’ 

-117.678 

104.5 

b t 

40617.6 

19662.8 








d t 

-14468.60 

-203.73 

-19.96 

2299.11 

.5885 

2632.15 

-20.069 

-43.53 



TABLE 3 Summary of Simulation Results 


Simulation 

Test 

Weighting 
Matrices 
fi. R 

Number of 
Samples, 

Regulator 
Performance 
Index, J 

u 

Samp lino 
Efficiency, ij 

Case 1 

1C, 

1 

1928 

.03315 

.1650 

Case 1 

10, 

1 

678 

.00337 

.3360 

Case 3 

50, 

1 

621 

.0256 

.0626 

Case 4 

10, 

1 

900 

.00279 

.3960 




Fig. 1. Typical output trajectories for a Gaussian Disturbance. 
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Fig. 2. (k) vs. k for 
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Fig. 3. The test input. 
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